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Abstract 

We look for the existence of stable (by strong interactions) tetraquarks, 
mesons containing two heavy antiquarks and two light quarks, specifically 
a \bbqq) state. This state is viewed as two \bq) particles bound by a poten- 
tial due to the exchange of the four light pseudoscalar and vector mesons, 
namely the 7r, n, p and uj. Couplings of these particles are obtained by the 
use of heavy quark symmetry. We find that the 1 = 0, J = 1 state is bound 
for a large range of parameters, while all other states are not bound (the 
1=1, J = system may be bound for some parameters at the extreme 
range of their allowed values). We discuss the reasons why the techniques 
used cannot be extended to mesons containing two c quarks or a c and b 
quark. 
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I. INTRODUCTION 



In addition to the usual three quark baryons and quark- ant iquark mesons, hadrons con- 
taining one or more heavy quarks may occur in other configurations, stable against de- 
cay by strong interactions. Baryons with one heavy antiquark and four light quarks, i.e. 
pentaquarks [|I], as well as baryons with two heavy antiquarks and five light quarks, i.e. 
heptaquarks 0, have been considered. Mesons containing two heavy antiquarks and two 
light quarks, tetraquarks have also been studied. The approach used in the last two 

references was to consider a tetraquark as two mesons, each containing a heavy antiquark, 
bound to each other through the exchange a 7r meson. 

In this work we shall look for bound states of mesons containing two b quarks and 
two light u or d quarks, namely a (bbqq) system. As mentioned above we shall view such a 
possible system as a bound state of two B mesons, two B* mesons, or a BB* combination. It 
turns out that the bound states are linear combinations of these configurations. The binding 
potential will be due to the exchange, in addition to the ir, of the 77, p and u> mesons. The 
relevant coupling constants are obtained in Section using heavy quark symmetry [|7| . The 
Schrodinger equations and their solutions are discussed in Section [TV|. We find that bound 
states may exists for the light quarks in a symmetric color combination (presented in Section 
PTT] ); this is the opposite of the conclusions from bag model arguments For a wide range 
of parameters we find that a bound state exists for the system with total isospin 1 = and 
total angular momentum J = 1. The range of binding energies below a possible decay to 
B and B* mesons is about 650 MeV to 250 MeV. The 1=1 and J = is unlikely to be 
bound. 

Binding of states containing two c quarks or a c and a b quark are discussed in Section |V]. 
Our methods cannot be extended to the study of these states. As Mo* — Mo ~ and for 
some charge combinations Mo* — Mo > m n we cannot use the multi-channel approach in 
which we neglect the mass difference of the heavy pseudoscalar and vector mesons. Diagrams 
where the ir is exchanged between a D and a D* result either in a weak very long range 
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potential or in an imaginary one, indicating that one should treat the DD* state as a three 
body DDn system. 



II. INTERACTION POTENTIAL 

We shall be interested in obtaining the potential between two heavy mesons. In the 
heavy quark limit |7j the pseudoscalar and vector mesons are degenerate in mass and the 
bound state problem has to be treated as a multichannel one. As in this limit the masses 
of the incident and final pairs are equal and the energy of each of these states is just the 
sum of the two heavy masses, we are also allowed to use non-relativistic methods. These 
approximations are adequate for mesons containing b quarks but, as we will discuss, are 
likely to break down for systems with lighter quarks. 

A. Heavy-light coupling constants 

As in Ref. || we shall use the formalism of heavy quark symmetry and of the nonlinear 
chiral model || to obtain the coupling constants for the exchange of light mesons. With B a 
and B* a denoting the fields for the pseudoscalar B and vector B*, a being an isotopic spin 
index, we can define a heavy, fixed four-velocity field H a , 

H a = ( 1+ ^" 7 ) {il*B a + 7 • B*;), H = 7 °tf V . (1) 

The free Lagrangian for H is 

C free = -iM VTr (Hd^H) . (2) 

The light Goldstone pions and the rj are described by a U(2) matrix 

X = exp j^(r-7r + ?7)j , (3) 

with F n = 132 MeV. Chirally invariant Lagrangians can be written with the help of the 
square root of the S field, £ = S 1//2 . We define two combinations 
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p, = ~ me - , « M = £ ce^ + e^o • (4) 

The interaction of heavy mesons with the light ones is given by 

C d = MdTr(H-f'y%H) . (5) 

The coupling constant d is not very well known; from a single pole fit to the data on D K 
and D — > K* semileptonic decays ||, one obtains d ~ 0.53. A straight forward use of the 
static SU(6), or more precisely SU{4) symmetry gives d = 1; reducing this by the ratio 
(tju/<7y)/(5/6) results in d = 0.7. The value d = 0.53, which we shall use, is a conservative 
estimate. If we include the strange quark content of the rj and rj — rf mixing the heavy 
meson- rj coupling has to be reduced to d'\ using the known 77 — rf mixing angle of 11° we 
find that d'/d ~ 0.4. 

The light vector meson fields, and u^, are combined as 

Pm = 72^' T + ^' ^ 
and we also define a field strength tensor 

F„ u (p) = d^p u - d v p^ - iglpp, p v \ . (7) 

The width of the p sets g = 3.93. The coupling of these vector mesons is accomplished by 
modifying £ free , Eq. (@) to 



Akinetic = ~iM VfjTr 



H{d lx - iagp^ - z(l - a)v^}H , (8) 



where is defined in Eq. (|J). This term is invariant under chiral transformations. Vector 
dominance suggests that a ~ 1, which is the value we use in this work. In addition we add 
a magnetic moment term 

-icM , 



C-c = Tr 

m v 



HYYF^(p)H\ . (9) 
m v is a fiducial light vector meson mass, which appears from dimensional considerations. 



The coupling constant c is constrained to [10 



2.27 < c/d< 3.75. (10) 
Phenomenology of the nucleon-nucleon potential fill requires the exchange of a + (0 ++ ) 



particle with a mass of ~ 550 MeV. This represent the effective potential due to the exchange 
of two pions in box diagrams with isobars in the direct channel and the exchange of the broad 
7r — 7r resonance /o(980). In our case, by considering both the scalar and vector mesons in a 
multichannel formalism, the box diagrams are automatically included and the exchange of 
the /o(980) contributes to that part of the potential whose range is shorter than that which 
we consider explicitly; hence we will not include any + (0 ++ ) exchanges. 



B. Interaction potential 

Using the coupling constants determined in the previous subsection it is straightforward 
to find the interaction potential. Let S* denote the spin of the light quark in the heavy 
meson, i = 1,2, while S and I denote the total spin and isospin of the two light quarks. The 
position space potential is a combination of central, spin-spin and tensor terms, 



V(x) = V Q {x) + (S 1 • S 2 )U 5 (x) + ( (S 1 • x) (S 2 • x) - ^S 1 • S 2 ) V T {x) , (11) 



where 
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(13) 



The 1/r 3 part of the tensor force, Vr, is too singular to permit a solution of the resulting 
Schrodinger equation; in addition there is a delta function at the origin lurking in that 
potential. As we have not taken into account the exchange of mesons whose mass is greater 
than that of the light vector mesons we will, as in Ref. ||, flatten the potentials at some r c ; 
we will choose r c = l/m p ; we will discuss how the results depend on r c for r c > l/m p . 

III. BOUND STATE CONFIGURATIONS 

In order to obtain a bound state we require that the long range part of the potential (tt 
exchange) be attractive. From Eqs. (|TT| , [12[) we see that the light quark configurations must 
be either I = l,S = 0orI = 0, S = 1; as expected the greatest attraction is in an orbital 
S'-wave state. As the space-spin-isospin of the light quarks is antisymmetric they must be 
in a symmetric color state, the 6; this implies that the two b quarks are in the 6*, forcing 
them to be in a spin singlet state. Thus the possible bound states are either a I = 1, J = 
state or a / = 0, J = 1 state. These are linear combinations of B's and B* y s: 

|/ = l,/ 3 ;J = 0) = i(l,/ 3 |l,m;l,m / } 



|/ = 0;J=l,J 3 = l) = ^(0,0|l,m;l,m / ) 

(\ B *m B *rn'} + \B* m B*™') - \B*™B m ') + \B m B*™')) , (14) 

with similar expressions for the other two J = 1 states. 

The expectation value of the mass operator in these states is twice the spin averaged 
mass 

3Mr. + Mr 

(M) = 2M eS = 2 6 B 4 + ■ (15) 
Eqs. ( |i~4l) tell us into what combinations of B's and _B*'s these states may decay. 

IV. SCHRODINGER EQUATIONS AND SOLUTIONS 
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A. Eigenvalue problem 



To order 1/M the Hamiltonian of this system is 

1 



H = 2M„ 



V 2 + V(x) . 



(16) 



A basis for the system is provided by the states I//3) \SSs) \lm), where l,m are the orbital 
angular momentum quantum numbers of the bound system. Note that l,m includes the 
orbital angular momentum of the heavy quarks as well. Since the potential is rotationally 
invariant, the total angular momentum 



J = S + L 



(17) 



is a good quantum number and can be used to label the eigenstates. The J = (singlet) 
state is written as 



(18) 



The J = 1 (triplet) state is a linear combination of \S — 1, I — 0) and \S — 1, I — 2) and 
can be written as 



^J=l- r ^1;0,1+ r ^1;2,1 > 



(19) 



where 3^j 3 / 5 is a generalized spherical harmonic obtained by combining spin and orbital 
angular momentum. Both the singlet and triplet states must further be supplemented by 
an isospin label, with I = or 1. 

With e = E — 2M e g, the eigenvalue equation for the singlet case is 

3, 



r = M eS (v --v s -6^ V; 



(20) 



prime denotes differentiation with respect to r. For the spin triplet case the tensor force 
couples I — to I — 2 and we get two coupled differential equations, 



u" = M. 



eff 



V + \v s -e)u + ^=V TW 



w" = M t 



eff 



M eff r^ 6 



_ _ y T — e j w _| — Vt« 



3^ 



(21) 
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B. Numerical results 



We now consider the numerical solutions to Eqs. (^H) and (|2ip for the energy eigenvalue 
e; d is fixed at 0.53 and the B, B* masses are taken to be 5.279 GeV and 5.325 GeV 
respectively, c/d is varied between 2.27 and 3.7. 

For the singlet case, we find that there are bound states only in the 1=1 channel, 
forcing the color wavefunction to be in the symmetric 6 representation, as mentioned earlier. 



We observe from the decomposition Eq. (fH|) of the singlet wavefunction that the lightest 



constituents the state can decay into is two B mesons. Hence the true binding energy is 

e B = 2M eff - 2M - |e| 

= 69MeV-|e|; (22) 

the state is bound if is negative. Binding occurs only for c/d> 3.6 and e# = 40 GeV at the 
extreme of the allowed range, c/d = 3.7. Hence we conclude that the singlet configuration 
is unlikely to be bound. 

We now move to the triplet configuration. Here, the lightest constituents the state can 
decay into is a B and a B* meson, hence the true binding energy is given by 

e B = 2M eff -M-M* - \e\ 

= 23 MeV - |e| . (23) 

Fig. 1 summarizes the numerical results for the allowed range of c/d. It is clear from the 
figure that there is a stable bound state over the entire allowed range of the coupling constant 
c. 

One interesting issue to explore is the sensitivity of the binding energy to the cutoff 
r c , below which the potential is flattened. Fig. 2 shows the dependence of on r c , for 
a central value of the ratio c/d = 3. The lightest particle above the vector mesons which 
can contribute to the potential is the /o(980), or the rj' (958). These correspond to a value 
of r c in the neighborhood of about 1 GeV" 1 . We see from Fig. 2 that a cutoff less than 

8 



about 1 GeV -1 leads to a very high binding energy; indeed, the non-relativistic formulation 
of the system breaks down. Hence the potential derived here cannot be trusted for values 
of the cutoff much less than 1 GeV -1 . We note that r c = l/m p is a conservative estimate 
for determining whether binding does or does not occur. 



We repeated the calculations for smaller reduced masses to see what happens for the 
D — D and D — B systems. From Fig. 3 we could claim that both of these will be bound 
for a triplet light quark spin configuration. This interpretation is too naive. As has been 
emphasized in Ref. [§] the heavy quark limit cannot be use for the study of binding of 
charmed mesons. For 7r exchange the range of the potential is nor l/m n , but with 



/i = Jm\ — (Mq — Md) 2 . This value is close to zero, and for several charge combinations 
negative. (For the B — B system this gives a 6% correction to the range of the 7r exchange 
potential.) Thus the breaking of heavy quark symmetry and even the breaking of isospin 
invariance will be crucial to the study of this system. As in the nucleon-nucleon case we 
could look at the single channel D — D problem. As there is no one-7r exchange potential for 
this configuration, this system is unlikely to be bound. The D — D* mass difference would 
force us to treat this as a three body D — D — n system. 
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FIGURES 

FIG. 1. Binding energy, eb, of the 1 = 0, J = 1 state as a function of the magnetic vector 
meson coupling c/d. 

FIG. 2. Binding energy, cb, of the 1 = 0, J = 1 state as a function of the cutoff for c/d = 3.03. 

FIG. 3. Eigenvalue of the Schrddinger Equation, e, of the 1 = 0, J = 1 state as a function of 
the reduced mass for c/d = 3.03 and r c = l/m p . 
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